Frictionless quantum quenches in ultracold gases: a quantum dynamical microscope 
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In this article, a method is proposed to spatially scale up a trapped ultracold gas while conserv- 
ing the quantum correlations of the initial many-body state. For systems supporting self-similar 
dynamics, this is achieved by implementing a many-body finite-time frictionless quantum quench of 
the harmonic trap which acts as a quantum dynamical microscope. 
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Research in ultracold strongly correlated states of mat- 
ter has recently been spurred by the experimental realiza- 
tion of quantum gas microscopes, allowing to detect with 
single site resolution and nearly unit-efficiency individ- 
ual atoms within a macroscopic sample in the strongly- 
interacting regime optical lattice [ , ]. These experi- 
ments are based on high-resolution optical imaging. In 
a complementary way, the non-equilibrium dynamics fol- 
lowing a quench of an external control parameter is of- 
ten exploited to probe quantum correlations in many- 
body systems [3]. Here we propose a scheme to imple- 
ment a quantum dynamical microscope, an engineered 
controlled expansion that allows to scale up an initial 
many-body state of an ultracold gas by a desired fac- 
tor while preserving the quantum correlations of the ini- 
tial state. The scheme rests on the possibility of driv- 
ing a self-similar dynamics in certain systems, which is 
a powerful tool to understand the evolution of quantum 
correlations. Scaling laws can often be exploited to de- 
scribe harmonically trapped ultracold gases, such as the 
Calogero-Sutherland model [ ], the Tonks- Girardeau [ - 
7] and certain Lieb-Liniger states [8], Bose-Einstein con- 
densates (BEG) [5, 9, 10], including dipolar interactions 
[11], strongly interacting ultracold gas mixtures [x..], and 
more general many-body quantum systems [13]. More- 
over, whenever the dynamics is not self-similar per se, it 
can often be assisted by tuning the interactions, either by 
means of Feschbach or confinement induced resonances, 
or time-modulation of the transverse confinement in ef- 
fectively low-dimensional Bose-Einstein condensates [ ] . 
Nonetheless, in spite of the scaling laws, the expansion 
dynamics in these systems generally induces undamped 
breathing of the cloud and distorts the quantum correla- 
tions of the initial state [5-13]. The method proposed 
in this Letter suppresses these effects in a finite-time 
non-adiabatic expansion which acts as a lens to zoom- 
up the initial state of the system. This is achieved by 
carefully engineering the time-modulation of the trap- 
ping frequency to induce a frictionless dynamics free from 
adiabaticity constraints. 

Self- similar dynamics- Let us consider a /^-dimensional 
many-body system composed of N indistinguishable par- 
ticles confined in a time-dependent harmonic trap, de- 



scribed by a Hamiltonian 
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is the D dimen- 
e{t) is a dimension- 



sional Laplacian operator, and e 
less time-dependent coupling constant which reduces to 
the identity at t = 0. We shall focus on systems with an 
interaction potential satisfying the relation 



V(Ax) = A"V(x) 



(2) 



under scaling of the coordinates. An equilibrium state ^ 
of the system (1) at t = 0, follows a self-similar evolution 



*({Xi},t) 



where /i is the chemical potential, and r{t) = 
Jq dt' /l)^{t'), whenever the scaling factor h = b{t) is the 
solution of the Ermakov differential equation 



b^uj'\t)b = ujl/b^ 



(4) 



with uq = cj(0), satisfying the boundary conditions 
6(0) = 1 and 6(0) = 0. This dynamics further requires 
the coupling constant and the scaling factor to be related 
by e{t) = which leads to the following cases: a) 

a = 2, e{t) = 1, with no need for auxiliary tuning of the 
interactions, as it happens for example in a quasi- ID 
Bose-Einstein condensate (BEG) in the Thomas-Fermi 
limit or a 2D Bose gas with contact interactions [ , ] , 
which has recently been pointed out as an instance of 
a quantum anomaly [ ]. b) a 7^ 2, e{t) = 6"~^ which 
does require external tuning of the interactions as in 
ID and 3D BEG to assist the self-similar dynamics 
(note that in a fast expansion the role of interactions 
might be disregarded, as usually done in time-of-flight 
experiments, and then the self-similar dynamics comes 
for free). If the initial state ^ is not in equilibrium, it 
will then follow the evolution as if the trapping potential 
is kept constant in the scaled coordinates and picking 
the overall phase in Eq. (3). Scaling laws manifest in 
non-local correlations of the gas such as the one-body 
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reduced density matrix (OBRDM) given by ^i(x, y;t) = 
Njd^2"- c^XAr^*(x, X2, . . . , xat; t)<l>(y, X2, . . . , xat; t), 
whose time evolution under self-similar dynamics can be 
conveniently written as [ 13] 



^i(x,y;t) = r^gi 



iVpo) 



exp 



i b x^ - y^ 
b ujQ 2/q 
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where Iq = ^/h/muJo^ and its Fourier transform n(k, t) = 
/ dxdy e*^"^^~^^^i(x, y; t), the momentum distribution 



n(k,t) 
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as well as any higher-order correlation function, i.e. the 
n-body reduced density matrix ^^({x^}^^^; {xj}^^^; t) = 



AT! 
{N-n)\ 
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N 



exp 



=n+l 
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dx,$*({x,}iIi)$({xaLi,{xJf^,+i;t) 

^n({t}Li;{^}Li;0). 



The Newton cradle experiment has singled out the ID 
Bose gas in the strongly interacting limit as a paradig- 
matic example of an integrable system where relaxation 
to equilibrium is suppressed even when perturbed to 
breakdown integrability [ ]. We shall illustrate our re- 
sults with a ID cloud of ultracold bosons in the limit of 
hard-core contact interactions, a Tonks-Girardeau (TG) 
gas, confined in a harmonic trap of frequency cjq, with 
single-particle eigenstates (/)n(^), and n = 0,1,2,... 
The many-body ground state of this system can be el- 
egantly described using an auxiliary wavefunction, a 
normalized Slater determinant [ ], . . . , xat) = 

;^^det^^~^^Q^j^n(^0? describing a spin polarized Fermi 
gas in the ground-state of the trap. This wavefunc- 
tion already includes the hard-core condition as a re- 
sult of the Pauli exclusion principle encoded in the de- 
terminant structure. The bosonic symmetry can be 
enforced by applying the antisymmetric unit function. 



A{xi 



<j<k<N 



e{xk — where e{x) = 1 



(-1) if X > (< 0) and e"(0) = 0. The Bose-Fermi map- 
ping relating both dual systems reads ^tg{xi^ . . . , xat) = 
A{xi^ . . . , f Ar)^F(^i, . . . , Xat). This is a highly non-local 
mapping, but being involutive, it leaves invariant any 
local correlation function such as the density profile, 
i.e. those quantities derived from the probability den- 
sity are shared by both dual systems. The computa- 
tion of non-local correlations remains a non-trivial task, 
but elegant expressions are known after the work by 
Pezer and Buljan [ ]. A many-body state in the TG 
regime obeys the self-similar scaling law in Eq. (3) 
whenever 6 is a solution of Eq. (4) [ ]. The scal- 
ing is clearly exhibited by the density profile of a TG 
gas riTGix.t) =NjdX2 . . .dXN\^TG{x,X2, . . . ,XN;t)\'^ - 
N J dx2 . . . dxN I (a;, X2 , . . . , xn ; ^) P 
its width being governed by b. 
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FIG. 1: Adiabatic Limit. The scaling factor following a 
quench for increasing quen ch time r eventually approaches 
the adiabatic result 7 = y^ujo/ujf) — vTO. The inset shows 
how b{tf) vanishes for large values of r (t/ = 20t). 



Breathing and correlation dynamics- Consider the self- 
similar expansion of a many-body system after quench- 
ing the trapping potential, between an initial uq and 
final Uf frequency in a finite quench time r, to scale 
it up by a factor 7 = (cjo/cj/)^/^. Two main features 
of the dynamics are to be tailored to create a quantum 
dynamical microscope: the subsequent breathing of the 
density profile and the evolution of quantum correlations 
gn- For illustrative purposes, we shall focus on those as- 
sociated with the OBRDM. First, we note that for the 
sudden expansion in free space {uj{t > 0) = ujf = 0), 
b{t) = y^Y^^~uj^ and for t (^0^^ b{t) ~ cjq^, b = ujq. 
Using the method of the stationary phase, n(k, t) ^ 
|27rcc;o^o/^l^^i(^ok/o/^7 ^ok/o/^)^ so the asymptotic mo- 
mentum distribution is the scaled density profile of the 
initial state [^'^, ^"^J, which in the TG regime can be fur- 
ther related to the momentum distribution of the dual 
system, the ideal Fermi gas [6, 7, 21, 22]. This mapping 
between local and non-local correlations is expected as 
the density decreases during expansion at a finite rate 
b. Indeed, within the scheme of symplectic tomography, 
the evolution under quadratic Hamiltonians leads to a 
dynamical covering of correlations in phase-space [23]. 
In an isolated system, a sudden quench of the trapping 
frequency between two given finite values also induces 
undamped breathing of the density profile. In the TG 
regime, this was shown by Minguzzi and Gangardt [ ], 
who further illustrated the dynamics of the momentum 
distribution, periodically oscillating between that of the 
initial trapped state and the quasi-momentum distribu- 
tion. Nonetheless, these two effects are ubiquitous in the 
family of systems supporting a dynamical scaling gov- 
erned by the Ermakov equation. For a non-vanishing 
cj/ / 0, increasing finite values of the quench time lead 
to a gradual suppression of the subsequent breathing, 
and one expects to recover the adiabatic limit whenever 
uj{t)/uj'^{t) ^1. As a specific example let us consider 
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FIG. 2: Adiabatic Limit. Momentum distribution of a Tonks 
Girardeau gas following a quench of the trapping potential 
[see Eq. (7)] with N = 10, uj{0) = ujq and oj{t) = cjo/10 for 
increasing quench times r =: A) 0.1, B) 1, C) 10, D) 100 
in units of cj^^. The dashed line represents the distribution 
following the quench, while the solid line corresponds to the 
desired adiabatic limit (t/ = 500, ko = 1/ muoo / h) . The insets 
show the evolution of the scaling function b{t) along and after 
the quench. The adiabatic limit is approached for increasing 
values of r ^ lOOcJo leading to a suppression of the breathing 
dynamics of the cloud and scaling of quantum correlations. 



the functional dependence of the trapping frequency to 
be given by 



uj{t) = ujQ -\- {ujf — ujq) tanh 



(7) 



with a characteristic quench time r. The resulting b{t) 
can be obtained by solving numerically the Ermakov 
equation subjected to the initial conditions 6(0) = 1, 
6(0) = 0. The effect of an increasing value of r on the 
scaling factor b{t) is shown in Fig. 1. Quenches of the 
potential in a finite time r ~ uJq^ distort the momen- 
tum distribution and induce breathing of the cloud. As 
it turns out, these quenches do lead to non-zero finite 
values of b{tf ^ r) and the scaling factor b{tf) deviates 
from the desired target value 7, as shown in Fig. 2. As 
r is increased and the adiabatic limit is approached, the 
variation of the scaling factor becomes slower and slower, 
and b{tf) 0. Eventually, the quench follows the adi- 
abatic trajectory b{t) = y^uj{t)^ for which the OBRDM 
reads 



^i(x,y;t) 



1 / X y 

b^^^ \b(t)'b(t)' 



(8) 



according to 

n(k,t) = b^{t) y^dxdy ^i(x,y;0)exp[z7k- (x-y)] 

= 6^(t)n(6(t)k,0), (9) 

and coincides with the initial momentum distribution up 
to the scaling factor b{t). Note that these expressions 
can be applied both for expansions {b{t) > 1) as well 
as compressions {b{t) < 1). Nonetheless, the required 
adiabatic time scale can be exceedingly long and we next 
tackle the problem of achieving a final scaled state in a 
predetermined time of expansion r. 

Frictionless dynamics: preserving quantum correla- 
tions during expansion- Following the theoretical propos- 
als in [15, 24] shortcuts to adiabaticity have been im- 
plemented in the laboratory both for thermal gases and 
Bose-Einstein condensates [^o, zu]. In the following, we 
show how shortcuts to adiabaticity can be exploited to 
control the dynamics of quantum correlations in the fam- 
ily of many-body systems given by Eq. (1). First, we no- 
tice that the Ermakov equation can be inverted to design 
a many-body fast frictionless trajectory between and ini- 
tial and a final trap in a given quench time r (a time ana- 
logue of the focal plane in an optical microscope), hence, 
providing a shortcut to adiabaticity. To this aim we en- 
force the scaling law to reduce to the initial state ^ at 
t = and its scaled- up form at t = r. This leads to the 
following boundary conditions for the scaling function, 
6(0) = 1, 6(0) = 0, 6(0) = 0, 6(r) = 7 = [cjqM]^/^ being 
the scaling factor, and 6(r) = 0, 6(r) = 0. This set of 
conditions can be used to fully determine the coefficients 
in an ansatz, say, of polynomial type b{t) = Xlj=o%^"^- 
One finds that 

b{t) = 6 (7 - 1) - 15 (7 - 1) + 10 (7 - 1) 8^ + 1, 

(10) 

with 5 = t/r, is a solution of the Ermakov equation that 
drives the evolution from the initial trapped state <^>(x, 0) 
to the final state ^(x, t) = 7"^/^^(x/7, t = 0) in a finite 
time r mimicking the adiabatic evolution. The required 
modulation of the trapping frequency can be obtained as 
well from Eqs. (4) and (10), and might involve imaginary 
frequencies associated with a repulsive potential cj^ < 
[ ] whenever the demanded expansion time is small, 
r ^ uJq^. Should that be required, the trajectory can be 
implemented in the laboratory as in [ ]. The upshot of 
the frictionless dynamics is that quantum correlations at 
the end of the quench (t = r, and only then) are those of 
the initial state scaled by a factor 6(r) = 7. In particular. 



^i(x,y;r) = -^gi ( -, -;0 
n(k,r) = 7^n(7k,0). 



X y 



(11) 



this is, it is the desired result of scaling the initial 
OBRDM. Similarly, the momentum distribution evolves 



Similar expressions hold for higher-order correlations, 
i.e. ^n({xi},{yi};r) = 7''^^n ({xi/7}, {yi/7}; 0), with 



4 



s = 0.00 
s = 0.25 
s = 0.50 
s = 0.75 
s= 1.00 




-30 -20 -10 10 20 30 

kb(t)/k^ 



In conclusion, for many-body systems supporting self- 
similar dynamics, we have shown how to scale up the 
system by means of a controlled expansion without mod- 
ifying the quantum correlations. As an alternative to the 
adiabatic dynamics, we propose to implement a fast fric- 
tionless quench of the trapping potential, which acts as a 
quantum dynamical microscope, leading to the scaled- up 
initial state at the end of the quench. 

Insightful discussions with J. G. Muga, M. B. Plenio 
and L. Santos are acknowledged, as well as financial sup- 
port by EPSRC and the European Commission (HIP). 



FIG. 3: (color online). Frictionless quantum quench of a 
Tonks- Girardeau gas. Evolution of the momentum distribu- 
tion exhibiting signs of dynamical fermionization at an inter- 
mediate stage of the expansion, before reaching the final time 
T at which it reduces to the scaled-up distribution of the ini- 
tial trapped state (7 = 10, N = 10). The inset shows the 
smooth evolution of the scaling factor in a frictionless expan- 
sion in an arbitrarily short quench time r, for both 7 = 
(solid line) and 10 (dashed line) expansion factors. 

{x^} = {x^}^^^. Moreover, as long as the initial state 
is an equilibrium state in the initial trap, so it is the 
state at r with respect to the final trap, preventing 
any non- trivial dynamics after the quench, for t > r 
if uj{t > t) = ujf. Nonetheless, at intermediate times 
t G [0, r) the momentum distribution exhibits a rich non- 
equilibrium dynamics, and can show for instance, evolu- 
tion towards the scaled density profile of the initial state. 
Note that this mapping is favored by a large scaling fac- 
tor 7. 

Implementing a frictionless dynamics allows to perform 
a controlled expansion towards the scaled-up state, as 
illustrated in Fig. 3. For a given final time r of ex- 
pansion, the momentum distribution approaches most 
closely the initial density around t ~ r/2. As shown 
in the inset, h has a maximum precisely at t = r/2 
around which it is approximately constant as in the 
asymptotic free expansion. For t > r/2 the evolution 
proceeds so as to reconstruct the initial state scaling it 
by the desired factor 7. We close noticing that the pro- 
cess is robust in the sense that within linear response, 
errors in the implementation of the trap modulation or 
a many-body perturbation represented by the operator 
£:W(x, t) (see the supplementary material) lead only to a 
quadratic decay of the fidelity between the target state 
<I>(r) and the resulting state at the end of the quench 
^^(r), T{r) = mrWir))\^ = 1 - ^^(r/rz)^ + 0{e^), 
where the Zeno-like time tz = ^/[(W^) - (W)^]^/^ with 
W = ^ Jo" Uo{-t)W{t)Uo{t)]dt where Uo is the time evo- 
lution operator. Hence, the effect of 6:W(x, t) is con- 
trolled by the ratio r/rz and can be suppressed in fast 
expansions. 
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SUPPLEMENTARY MATERIAL 

STABILITY AGAINST DEVIATIONS FROM THE 
DESIGNED TRAP MODULATION 

Using the Ermakov equation Eq. (4), and the scahng 
factor b{t) in Eq. (10) one can obtain the trajectory 



60(7 - l)s(s(2s - 3) + 1) ((7 - l)(3s(2s - 5) + 10)s^ + l) 



((7 - l)(3s(2s - 5) + 10)s3 + ifw^rZ 



(12) 



with s = t/r, which leads to a successful implemen- 
tation of the quantum dynamical microscope. Clearly, 
the pair {b{t)^uj{t)}) is not unique, and provided that 
b{t) satisfies the boundary conditions at t = 0,r, many 
alternatives can be found. In practice, deviations of 
these trajectories are likely to occur, and we next con- 
sider the stability against a time-dependent frequency 
shift S{t) associated with the one-body perturbation 
£W(x,t) = e^,v{xi,t) = sEi^rn5u^{t)xl {5uj\t) = 
2uj{t)S{t) + As an illustrative example, let us con- 

sider the fidelity J^{r) = |/(t)P of a N-body state in the 
TG regime which at t = is the ground state of the trap, 
and where /(r) = {^{r)\^'{r)) is the fidelity amplitude 
between the target ^(r) and perturbed ^\r) states. 

In an expansion on the strength of the perturbation, 
the Zeno-like time governing the decay of J^{t) = 1 — 
e^{r/rz)^ + 0{e^) is found to be {N > 1) 



rz 



N dtSuj^{t)b^t)' 



(13) 



According to Eq.(lO), b{t) < 7 for ah t e [0,r], so that 
a high-fidelity microscope requires 



r < 



2V2UJ0 



(14) 



Hence, the fidelity decay is governed by the time average 
of the square of the frequency shift 5uj'^ = ^ dtSuo'^ {t) . 
In particular, for a Gaussian white- noise stochastic vari- 
able 5{t) of zero mean and covariance Jq, averaging over 
different realizations leads to {5uj'^) = ^q. We further 
note that the robustness against shifts in time of the 
ideal trajectory co'(t), is warranted by the smoothness of 
the scaling factor b{t) in the final stage of the quench, in 
the sense that b{t) 7, b{t) :^ 0. 



STABILITY AGAINST PERTURBATIONS OF 
THE MANY-BODY HAMILTONIAN 

Physical implementations of the family of ultracold 
gases described by Eq. (1) generally include deviations 
represented by a many-body perturbation W(x). In the 
following, we study the robustness of the quantum dy- 
namical microscope along a trajectory b{t) designed ig- 
noring a pair-wise perturbation £W(x) = e'^-^j v{xi — 
Xj), and compare the result with that of the full dynam- 
ics governed by the perturbed Hamiltonian. To illustrate 
this we shall consider and ultracold cloud in the Tonks- 
Girardeau regime as described in the text, for which the 
trapped modulation uj{t) and scaling function b{t) are 
designed. Assume that the real cloud prepared in the 
laboratory is a strongly interacting Bose gas, but not 
strictly in the TG regime. Generally, a Lieb-Liniger Bose 
gas with finite interactions in a time-dependent trap, de- 
scribed by a many-body Hamiltonian 



N 

E 

i=l 



2m' 



1 



Kj 



does not follow a self-similar expansion under a trap 
modulation ijo{t) unless the interaction strength is tuned 
simultaneously according to gfj){t) = g^^/b{t) as dis- 
cussed in the main text, case b) with a = 1. The effect 
of finite interactions can be taking into account by consid- 
ering the perturbation on the TG dynamics represented 
by the two-body J'^-pseudopotential [ ], 

i<j 



The breakdown of self-similar dynamics induced by it can 
be suppressed by imposing 



r < 



2h?'kl Cn 



(16) 



where Cn < N'^ is a coefficient depending on the number 
of particles. 



The upshot is that perturbations render unfeasible the quench time r associated with the expansion, 
an adiabatic version of the quantum dynamical micro- 
scope, while employing an engineered finite- time expan- 

sion makes it stable against both one and two-body per- 
turbations. This follows from linear response theory 

which to second order on the strength of the perturbation ^rsiT^Jaoir'''^' ^ '''' '^^^ ^' 

imposes a quadratic decay of the fidelity controllable by ' 



